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ABSTRACT 

The local luminosity function at 25 /zm provides the basis for interpreting 
the results of deep mid-infrared surveys planned or in progress with space 
astrophysics missions including ISO, WIRE and SIRTF. We have selected 
a sample of 1458 galaxies from the IRAS Faint Source Survey with a flux 
density limit of 250 mJy at 25 fim. The local luminosity function is derived 
using both parametric and non-parametric maximum-likelihood techniques, and 
the classical l/V max estimator. Comparison of these results shows that the 

estimate of the luminosity function is significantly affected by the Local 
Supercluster. 

A maximum-likelihood fit to the radial density shows no systematic increase 
that would be caused by density evolution of the galaxy population. The density 
fit is used to correct the l/V max estimate. We also demonstrate the high quality 
and completeness of our sample by a variety of methods. 

The luminosity function derived from this sample is compared to previously 
published estimates, showing the prior estimates to have been strongly affected 
by the Local Supercluster. Our new luminosity function leads to lower estimates 
of mid-infrared backgrounds and number counts. 

Subject headings: infrared: sources - luminosity function 
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1. Introduction 

Much of the effort to study infrared-luminous galaxies has centered on wavelengths 
greater than 50 fim. Modeling work is focused on the near-IR (e.g. Chokshi et al. 1994) 
and far-IR (e.g. Hacking et al. 1987; Rowan-Robinson et al. 1996) portions of the galaxian 
spectrum. However, the mid-infrared is well-suited for studying starburst and ultraluminous 
galaxies. About 40% of the luminosity from starburst galaxies is radiated from 8-40 /im 
(Soifer et al. 1987). Extinction effects are small, and problems due to infrared cirrus are 
minimized. Most importantly for space astrophysics, for a fixed telescope aperture, the 
spatial resolution is higher at shorter wavelengths, and the confusion limit lies at higher 
redshifts. Recent work using the Infrared Space Observatory (ISO) (e.g. Knapp et al. 1996; 
Boulade et al. 1996; Rowan-Robinson et al. 1996) shows the relative importance of the 7 
/zm and 15 /im bands for galaxy studies. The Wide-Field Infrared Explorer (WIRE), a Small 
Explorer mission due to launch in late 1998 (Hacking et al. 1996; Schember et al. 1996), 
will conduct a very deep survey at 24 /im to study starburst galaxy evolution. The Space 
Infrared Telescope Facility (SIRTF) is also expected to conduct surveys in mid-infrared 
bands. To interpret the results of these surveys now in progress or soon to commence, it is 
necessary to better understand the mid-infrared properties of galaxies in the local Universe. 

The 25 /im luminosity function provides the basis for predicting the faint source counts 
in the mid-infrared. The empirical model of Hacking & Soifer (1991) uses an analytic fit 
to the luminosity function derived by Soifer & Neugebauer (1991). This function was 
estimated from a complete subsample of the Bright Galaxy Sample (Soifer et al. 1987) 
containing 135 galaxies to a flux density limit of 1.26 Jy. The availability of many more 
redshifts of IRAS galaxies (principally from the 1.2 Jy Survey (Strauss et al. 1990; Fisher 
et al. 1995)) enables a much larger sample to be studied, reducing uncertainties at high and 
low luminosities. 
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In this paper we present the selection of a large galaxy sample that is flux-limited 
at 25 /im, and derive the local luminosity function based on this sample. The sample 
selection is described in the next section. In Section | we describe the l/V max and the 
maximum-likelihood estimators for deriving the local luminosity function, and present the 
results. The completeness of the sample is discussed in Section £|. In Section |5] we calculate 
the radial density distribution of the sample using a maximum-likelihood method. The 
radial density fit is used to correct the l/V max estimate of the local luminosity function, 
as well as the redshift distribution with which the luminosity function can be compared. 
Section 6 includes discussions of the different luminosity function estimators, a comparison 
of our newly derived luminosity function with previous estimates, the implications for 
mid-infrared backgrounds and number counts, and the effects of evolution on the derivation 
of the luminosity function. The color properties of the sample are treated in another paper 
(Fang et al. 1998). 

2. Sample Selection 

Our sample is based on a selection from the IRAS Faint Source Survey (FSS; Moshir 
et al. 1992). The principal data product of the FSS is the Faint Source Catalog (FSC). 
The FSC was produced by coadding IRAS detector scans before extracting sources, and is 
roughly one magnitude deeper than the Point Source Catalog. Another FSS database is the 
Faint Source Reject File (FSCR), which contains possible detections that were not included 
in the FSC for assorted quality-control reasons. This section treats the selection of sources 
from the FSC and the FSCR in turn. 

To minimize the effects of cirrus and contamination from stars, we chose a sky coverage 
limit of > 30°. In addition, the same excluded areas used by Strauss et al. (1990) were 
excluded from our sky coverage region. At high Galactic latitudes, the excluded regions 
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consist almost wholly of gaps not covered by IRAS, and the Large and Small Magellanic 
Clouds. Our sample covers 47.7% of the sky. 

The 90% completeness limit for the FSC lies at a 25 /mi flux density of 210 mJy for 
sky covered by 2 HCONs, and at 170 mJy for coverage of 3 HCONs, for |6| > 10° (Moshir 
et al. 1992). Our preliminary searches of redshift databases indicated, however, that for 
such faint flux density limits, many of the faintest objects would have no redshift available. 
A flux density limit of 250 mJy was found to have a small redshift incompleteness (see 
below) and was selected for our sample. Furthermore, sources with moderate (SNR=3 to 
5) or high-quality (SNR > 5) flux densities at 25 /mi are included in the sample. 

A set of weak color criteria were chosen to further discriminate against stars. The great 
majority of galaxies have larger flux density at 60 /mi than at 25 /mi. For these sources, 
no constraint was placed on the 12 /mi - 25 /im color. A small percentage of galaxies will 
have a smaller flux density at 60 /mi than at 25 /im. For these sources, we have limited the 
ratio of F u (25) / F V (6Q) to lie between 1 and 1.6, and have further constrained F„(12) / F v (25i) 
to be less than 1 to avoid admitting large numbers of stars into the sample. Color-color 
diagrams of the sources in the sample that have been identified as galaxies with redshifts 
(see below) are displayed in Figure [l| These diagrams show that only a handful of galaxies 
satisfy F u (25) / F u (60) > 1, and that our color criteria are fairly liberal for sources identified 
as galaxies. 

To summarize, our sample selection criteria are: 
F„(25) > 250mJy; 

moderate or good quality detection at 25/im; 
\b\ > 30°, and not in Strauss et al. excluded zone; 
F v (25)/F u (60) < 1 and no constraint on ^(12)7^(25), 
or 1 < ^(25)7^(60) < 1.6 and F v (12) / F v (25) < 1. 
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1619 sources in the FSC meet the above criteria. To identify galaxies and non-galaxies 
(stars, infrared cirrus, etc.), and to obtain redshifts for identified galaxies, a positional 
match of these 1619 sources was made with the 1.2 Jy Survey catalog, the November 1993 
public version of J. P. Huchra's ZCAT, the NED database, and the SIMBAD database. A 
matching radius of 1 arcminute was used to match the FSC sources with these catalogs. 
Additional redshifts were culled from J. P. Huchra's private ZCAT catalog, and from 
additional observations made at the Fred L. Whipple Observatory (FLWO) 1.5m. This 
matching process resulted in 1438 galaxies with redshifts. Of the remainder, 166 sources 
were identified as Galactic objects or HII regions in nearby galaxies, leaving 15 sources 
identified as galaxies without redshifts. 

The selection of sources from the FSCR is more complicated due to the low reliability 
of sources in that database. Applying our selection criteria to the Faint Source Reject File 
results in 851 sources. These were matched with the IRAS OPTical IDentification (OPTID) 
database at IPAC. The OPTID database is a special version of the FSC and FSCR with 
optical identifications added from the Guide Star Catalog (GSC: Lasker et al. 1990), the 
Tycho Input Catalog (TIC: Egret et al. 1992)), and the COSMOS/UKST Southern Sky 
Object Catalogue (COSCAT: Yentis et al. 1992). Sources matching with stars within la 
error ellipses on the OPTID Plots were excluded. We examined those sources which have 
no match with either a star or a galaxy within la error ellipses by studying the OPTID 
Reports and by comparing with the optical Digitized Sky Survey images. We found only a 
few of these sources to be possible galaxies which, along with the reject sources that match 
with galaxies within the la error ellipses, gives a small sample of 27 sources. Among these, 
there are 5 sources which are extremely blue, and have null SNRs at 60 /im. ADDSCANS 
(one-dimensional coadds of IRAS scans) of these sources show no obvious detections at their 
locations at all four IRAS bands. No redshifts for these sources were found when matching 
them with the redshift catalogs (see above), and we suspect that they are not real galaxies. 
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The final yield from the Reject File is 22 sources. 

Matching the 22 FSCR sources with the same redshift catalogs gives 20 sources with 
redshifts. To verify that these matches are not random, a sample of 851 sources randomly 
distributed in our sampling regions was generated. The number of matches with galaxies is 
only 2. This indicates that these 22 reject sources are not likely to be just random matches 
with galaxies. Furthermore, it shows that the number of coincidental matches in our total 
sample should also be very small. 

The breakdown of sources for redshifts for our entire sample (FSC plus FSCR sources) 
is listed in Table [l]. The previously unpublished redshifts obtained at FLWO are tabulated 
in Table |2|. The total number of galaxies with redshifts is 1458. There are 17 galaxies 
remaining without redshifts. 

The 1.2 Jy Survey catalog also provides the ADDSCAN flux densities for the sources 
extended at 60 /xm. For those sources which have no ADDSCAN flux densities in the 1.2 
Jy Survey but are extended at 25 ^m, we have obtained new ADDSCAN flux density 
measurements and substituted these for the FSC values. 

Figure |2] shows the log(iV) versus log(i^(25)) plot for the sample of galaxies with 
redshifts. The points follow a slope of -1.5 all the way to our flux limit, demonstrating the 
high completeness of our sample as a function of flux. The apparent excess at high flux 
values is most likely due to the Local Supercluster. Further checks of completeness are 
presented in Section |], including an investigation of possible systematic effects of flux errors 
on derived luminosity functions. 

For each galaxy, secondary distances are used when available. Otherwise, distances 
have been assigned using a linear Virgocentric inflow model (Aaronson et al. 1982) and a 
Hubble constant of 75 km s _1 Mpc -1 (see Appendix [A] for the case of a different Hubble 
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constant). The same criteria as Soifer et al. (1987) have been used to identify Virgo cluster 
galaxies (with a distance of 17.6 Mpc). The monochromatic luminosity of each source is 
computed from the distance, using a /^-correction that assumes a power-law slope of the 
SED between 12 /zm and 25 /xm. The monochromatic luminosities are expressed as vL v 
and have units of solar luminosities. 



3. Luminosity function results 

3.1. 1/V max and parametric maximum likelihood 

The classic method of estimating luminosity functions is via the 1/V max estimator 
(Schmidt 1968). This technique is non-parametric, so no analytic form of the luminosity 
function need be assumed, but the data must be binned. The infrared luminosity functions 
presented by Soifer et al. (1987) and Soifer & Neugebauer (1991) were derived via this 
method. The space density Pl and its error a p are computed from the following quantities: 



max 



^)(E^) 2 , m 

v ' \ max / 

where Q is the solid angle of the survey, V max is the maximum volume to which the object 
could have been detected, and the sum is over all galaxies in each luminosity bin. The 
assumption that galaxies are distributed uniformly in space may be tested by checking that 
V/Vmax is 0.5 in each luminosity bin. 

The luminosity function may also be derived via maximum likelihood techniques 
that are independent of density variations (e.g. Sandage, Tammann, & Yahil 1979). The 
parametric method has the advantage that no binning of the data is required, but an 
analytic form of the luminosity function must be assumed. For all of the parametric 
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maximum likelihood fits in this paper, we have assumed a cumulative luminosity function 
of the form used by Yahil et al. (1991) in the analysis of the 1.2 Jy Survey: 



with its corresponding differential luminosity function 

* (i) = (i + i^x)* (i) - (4) 

A flux-limited sample contains only a small number of very sub-luminous galaxies, 
making it difficult to determine the luminosity function at these luminosities. We have 
followed Yahil et al. (1991) in imposing a lower limit on luminosity L s = 47rr^z// m , where r s 
is the distance corresponding to a velocity of 500 km s _1 , and f m is the survey flux limit. 
For f m = 250 mJy, L s is 4.0 x 1O 7 L . Then the minimum detectable luminosity at distance 
r is L min (r) = Max(L s , 47ir 2 vf m ). The probability of detecting a galaxy of luminosity L; t at 
distance r, ; is then 

min , . 

(5) 

otherwise 



The likelihood function is the product of all of these probabilities: 

A=-2j>/(L i |r i ). (6) 

i 

This (reciprocal) likelihood function is minimized to determine a, (3, and L*. The variance 
of these parameters is estimated by making use of the asymptotic normality of the maximum 
likelihood estimators for a large sample (Kendall, Stuart & Ord 1987). 

The likelihood function contains ratios of the differential and cumulative luminosity 
functions, so another calculation is required to find the normalization C. One method is 
to use the n\ estimator of Davis & Huchra (1982) after the shape parameters have been 
determined. The quantity n\ is the number density of galaxies with L > L s . (See Equations 
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ll and 13 of Yahil et al. (1991) for the formulas appropriate to our choice of the parametric 
luminosity function.) An upper limit to the redshift must be chosen to fix the sample 
volume used in this density estimator. Based on a calculation of n\ in radial shells, we have 
chosen an upper limit on redshift of 20,000 km s -1 (z ~ 0.067) for the normalization of the 
maximum likelihood fit. (However, this upper redshift limit is not applied in determination 
of the shape of the luminosity function.) Effects of incompleteness or evolution may become 
significant beyond this redshift. A maximum-likelihood fit to the density variations of 
the sample as a function of redshift (see Section further indicates that the density is 
well-behaved to this limit. A second means of determining the normalization is to require 
that the number of sources predicted by the luminosity function over some redshift range 
equals the number in the sample. Computing the normalization from source counts in this 
way from 500 km s _1 to 20,000 km s _1 gives essentially the same normalization as the n\ 
estimator. 

The derived parameters for the maximum likelihood fit are listed in Table |^. Luminosity 
functions derived by the 1/V max method and the parametric maximum likelihood method 
are shown in Figure |3[ We show the luminosity functions in the form of the visibility 
= $L 2 5 . The visibility is proportional to the number of galaxies visible per log luminosity 
bin, and it has the advantage of showing the relative numbers of galaxies in a flux-limited 
sample. The open triangles are from the 1/V max estimate of the entire sample, and the 
solid curve is the parametric maximum-likelihood fit to galaxies with redshifts greater 
than 500 km s _1 , with normalization from a volume-limited sample as described above. 
The solid curve is truncated at L s = 4.0 x 10 7 L Q . The 1/V max estimate (made in bins of 
\og{vL u ) = 0.4) was converted to visibility by multiplying p L by L L5 j|^j. 

A plot of V/V max as a function of luminosity for our sample is included in Figure |j. 
The variations in this statistic for the full sample and northern sample at luminosities less 
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than 10 10 L Q are most likely due to the Local Supercluster. The southern sample is much 
more uniform at these luminosities, supporting this interpretation. 

The general agreement between the 1/V max luminosity function and the shape of the 
parametric fit in Figure 0, together with the uniformity of the V/V max statistic, confirm 
that the parametric fit is a good description of the shape of the local luminosity function. 
The 1/V max point at \ogvL v = 8.6 corresponds to the luminosity at which the flux limit 
samples the Virgo cluster, and hence lies well above the solid curve. 

The number of galaxies without redshifts (17) is slightly more than one percent of 
our sample. One method for estimating the effects of this known incompleteness on the 
luminosity function is to assign to each of these galaxies the median redshift of their log flux 
bin, and recalculate the luminosity function. When the median redshifts are computed for 
the whole sample, all the assigned redshifts are around 7,000 km sec -1 . When the median 
redshifts are calculated from those galaxies whose redshifts were not obtained from public 
sources (a more realistic estimate), the assigned redshifts range from 12,000 to 18,000 km 
sec -1 . Since all these galaxies have a 25 /im flux density within a factor 1.5 of the flux limit, 
their luminosities all fall near L*. When the luminosity function is recalculated and plotted 
as a visibility, we have found that the asymptotes are unchanged, and the peak visibility is 
increased by a few percent. 



3.2. Nonparametric maximum likelihood 

A non-parametric maximum-likelihood method (a.k.a "Stepwise Maximum Likelihood") 
has been used to calculate the galaxy luminosity function in various surveys (Lynden-Bell 
1971; Choloniewski 1986; Efstathiou et al. 1988; Saunders et al. 1990; Marzke et al. 1994). 
This method does not assume a specific functional form for the luminosity function, but 
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rather describes it by a non-negative step function (i.e. requires binning) which makes 
the observed galaxy luminosity distribution as likely as possible. We have followed the 
prescriptions set out by Efstathiou et al. (1988) and Saunders et al. (1990) for calculating 
the non-parametric luminosity function. 

As in the parametric maximum-likelihood method, the non-parametric approach also 
eliminates the effect of any density inhomogeneities, but loses the overall normalization 
by assuming that the galaxy spatial and luminosity distributions are independent. The 
normalization of the non-parametric function has been chosen to give the same source-counts 
as the parametric function over the same bins. 

Figure |5| shows the visibility function obtained by using the non-parametric maximum- 
likelihood method (circles), which agrees within the error bars with our parametric results 
(solid line, from Figure |3|). The non-parametric method is commonly used to assess whether 
the parametric functional form is acceptable. Figure |5] shows that the parametric form 
of Figure |^ gives a reasonable description of the density- independent 25 fim luminosity 
function. We have found that the non-parametric function is sensitive to the bin width, 
however. Figure [5] shows the case for bins in logL = 0.2, but for a bin width of 0.4, 
the non-parametric results is systematically higher than the parametric curve for high 
luminosities. The coarser binning probably results in a loss of information and a less 
reliable luminosity function. The choice of bin size is a compromise between this effect and 
tolerating the statistical uncertainty from having too few sources in each bin. 



4. Completeness Checks 



Systematic variations of noise across the sky may change the detection quality and 
cause incompleteness problems for the sample. In the 25 /^m passband of IRAS, noise 
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is dominated by photon shot noise from the zodiacal background, and hence is a strong 
function of the ecliptic latitude. We have examined the noise variation in the sample 
and in the corresponding FSS plates and found two regions that contain higher noise at 
25/mi. They are essentially the two quadrants defined by (0° < / < 180°, b < —30°) and 
(180° < / < 360°, b > 30°), which contain the ecliptic plane. Most of the moderate-quality 
detection sources (generally, those with SNRs between 3 and 5) in our sample are in these 
two regions, indicating the systematic effect. 

To estimate this systematic effect on our sample completeness, we compared luminosity 
functions estimated for two subsamples of our 250-mJy-limited full sample. The first 
subsample is defined by raising the flux-density limit to 400 mJy in the two quadrants 
of the sky (0° < / < 180°, b < -30°) and (180° < I < 360°, b > 30°) where the noise is 
systematically higher, and keeping the 250 mJy limit in the other two quadrants. Most 
of the moderate quality detection sources are left out of this subsample, so we call these 
1050 sources the high-quality subsample. The second subsample consists of those 408 
sources in the full sample that are left out of the high-quality subsample, and is therefore a 
low-quality subsample. Figure ^| compares the luminosity functions of the high-quality and 
the low-quality subsamples, using the l/V max estimator since it accomodates the different 
flux limits more easily than either of the maximum-likelihood estimators. There is no 
systematic deficiency of galaxies in the low-quality subsample. The luminosity functions 
agree within the la error bars, with the exception of the point at \ogvL v = 8.6. The 
low-quality subsample includes galaxies in the Virgo cluster with 25 /j,m flux densities 
between 250 and 400 mJy, most of which fall into this luminosity bin. 

In Figure [7| we compare the l/V max luminosity functions of the high-quality subsample 
and the full 250 mJy-limited sample. They agree very well, indicating that the full sample is 
complete at 250 mJy flux-density across the entire sampling sky. Therefore we are confident 
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that the full 250 mJy-limited sample is well-defined. 

As another check of the completeness of our full sample, we have computed the 
cumulative probability of finding a galaxy of luminosity L at distance r, following the 
method of Yahil et al. (1991). This probability is the integral of Equation |5] and is 
F(L\r) = ^/(L)/^/(L min (r)). This probability should be uniformly distributed on the 
interval [0,1] independent of distance. If the sample were incomplete at the faint end, we 
would expect a decrease in the relative number of galaxies as F — > 1. 

We have computed F(Li\ri) for each galaxy in the sample, then computed the histogram 
of this statistic in different distance ranges, normalized to 1 in each range. The result is 
shown in Figure || There is no systematic decrease in the relative galaxy counts as F — > 1 
(albeit the statistics are poor beyond 400 Mpc). This result provides further confirmation 
that our 250 mJy-limited sample is complete to these distances. Since this method relies 
upon ratios of the luminosity function, it should be insensitive to inhomogeneities in the 
galaxy distribution. 



5. Density Variations and the Corrections for V/V max 

Just as the galaxy luminosity function independent of the inhomogeneities of the 
density field can be obtained, the spatial density variations independent of the galaxy 
luminosity function can also be calculated using the maximum-likelihood method (Saunders 
et al. 1990). Here we only investigate the radial density distribution. The probability for 
detecting a source of luminosity Li at radial distance is 

= p{n)dVi 

Io max,t p{r){dV/dr)dr' K ' 

where p(r^) is the density at r«, and r max ^ is the maximum distance for the source to 
be included in the sample. The likelihood function for all the detected sources would 
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be C = Y\pi, given the density distribution p{r) we are trying to calculate. In the 
non-parametric approach, this distribution can be found by maximizing InC with respect 
to the binned density field p(rj), which gives 

p(^) = ^(£^r\ (8) 

where V e ff,i = X)I=cT '* p( r k)Vk, Nj is the number of sources in the j-th distance bin with 
volume Vj, is the fraction of the bin contained in [0,r moa; i] where the i-th source is 
observable, and the sum is over all sources. 

Figure ||] illustrates the non-parametric maximum-likelihood estimate of the radial 
density variation (solid line) as compared with the flux-limited observed distribution 
(dashed line). Here the radial distance bins are in logarithmic scale, and the density 
variations are normalized to the average density over all bins of interest. The method 
correctly recovers the density enhancement at ~ 17 Mpc caused by Virgo. 

We do not find a systematic increase in the radial density at large distances. If our 
sample is assumed to be complete, this result means that there is no indication of density 
evolution in this sample. We return to this point in section |67|. 



The benefits of calculating the radial density variations in our sample are two-fold. 
First, it may be used to correct the V/V max statistic and the traditional 1/V max estimate of 
the density-dependent luminosity function (Schmidt 1968; Felton 1977). For the V/V max 



test, since 



, r r p{r')D 2 dr' . . _„ , 

j o jf ax p(r>)D*dr>Pv> IJ dr /p ma3: (Jp pif) D 2 dr')d(J^ p(r)D 2 dr) _ 1 

^ max p{r)DHr ~ {§ max p(r)D 2 dr) 2 ~ 2' ^ 

where p(r) is the radial density field, r max is the maximum comoving distance for a source 
to be included in a sample, and D is the effective distance (Longair 1978), we can modify 
the V/V max test so that 

V _ r oP {r)D 2 dr 



V m ax Io max p(r)D 2 dr 
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which gives V/V max = 1/2 for an inhomogeneous radial distribution. This also suggests 
that we use the V e ff,i m equation § instead of V max in the l/V max estimator to calculate a 
density-corrected luminosity function. The p(r) calculated from the maximum-likelihood 
method can provide these corrections. The results are shown in Figure [TO for the new 



V/V max values and in Figure [TTJ for the 1/V max luminosity function. Compared with the 
results without the inhomogeneity correction in Figure £|, the V/V max values are closer 
to 0.5 for those bins most strongly affected by the Local Supercluster. (However, the 
corrected V/V max values are not suitable as a test of the sample completeness, since the 
density correction will mask any incompleteness in the sample.) The new 1/V max luminosity 
function is smoother and closer to the ones estimated by the maximum-likelihood method. 
As mentioned in Saunders et al. (1990), the binning of the distance in Figure § may account 
for the residual effects of the density inhomogeneities in these Figures. 

Second, the "true" density variations in Figure ^| also correct the observed redshift 
distribution. The density-variation-corrected redshift distribution should be predicted by a 
density-independent luminosity function estimated in the same sample. This provides an 
independent test to verify the luminosity function. Figure shows our results. Here the 
dashed histogram is the observed redshift distribution. (Here we show the source-counts 
in distance bins instead of redshift bins since some galaxies have secondary distances.) 
It is divided by the density of the corresponding distance bins in Figure |] which gives 
the solid-line histogram. The source-counts affected by large-scale structures, such as the 
ones near Virgo, are corrected by this procedure. The dotted line shows the prediction 
by the parametric density-independent luminosity function of Figure ^| and Table 0. This 
prediction fits the solid-line histogram well, demonstrating the reliability of our luminosity 
function. Both the solid-line histogram and the dotted line are normalized so that they 
contain the same total number of sources as observed within 20,000 km s _1 . 
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6. Discussion 

6.1. Limitations of the methods 

The preceding sections have presented the 25 /im luminosity function calculated by a 
number of methods. In this section, the limitations of the various methods are revisited in 
greater detail. 

The 1/V max estimate of the luminosity function varies the most from the other 
estimators, particularly at luminosities less than L#, where our sample is most strongly 
influenced by the Local Supercluster. As noted above, this estimator assumes that the 
3-dimensional galaxy distribution is homogeneous. The error bars in all of our l/V max 
estimates similarly depend on the assumption of uniformity, and do not reflect the additional 
uncertainties that arise from an inhomogeneous distribution. We consider the 1/V max 
estimates to be the least reliable descriptions of the 25 /im luminosity function out of those 
we have presented. 

The parametric and non-parametric maximum-likelihood estimates should be 
independent of structures in the galaxy distribution. However, these estimators assume 
that the luminosity function is independent of density. This assumption is difficult to 
check. In the case of a large 60 /zm-selected sample, Yahil et al. (1991) demonstrated that 
their similarly-derived luminosity function provided a good description of the luminosity 
distribution of galaxies even in regions of high density. 

The non-parametric maximum likelihood method applied to our sample is sensitive to 
the choice of binning, as described in section |3.2| . We therefore consider the parametric 
maximum-likelihood description of the 25 fim luminosity function with parameters in Table 
H] to be the best out of those we present in this paper. 
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6.2. Comparison with other luminosity functions 

As mentioned in the Introduction, Soifer & Neugebauer (1991) derived the local 
luminosity function at 25 fim from a complete subsample of the 60-/im-selected Bright 
Galaxy Sample (Soifer et al. 1987). This subsample contained 135 galaxies to a flux density 
limit of 1.26 Jy, and the luminosity function was determined using the classical 1/V max 
estimator. Their luminosity function is plotted as a space density in Figure [13], together 
with our parametric and 1/V max estimates. 

Our flux density limit of 250 mJy samples the Virgo cluster at \ogvL u = 8.6 (section 



3T| ), and our 1/V max estimate is much higher than the parametric curve in that bin. The 
Soifer & Neugebauer points in Figure |3| in log vL v = 9.4 bin is also much higher than 
the curve, as expected since their flux density limit is five times higher than ours. We 
conclude that the Soifer & Neugebauer luminosity function is strongly affected by the Local 
Supercluster, particularly because the Bright Galaxy Sample was selected to be visible from 
Palomar Observatory (Soifer et al. 1984) and does not include the south Galactic cap. 



6.3. Implications 

The improved luminosity function that is the main result of this paper has implications 
the luminosity density of the local Universe, infrared number counts predictions at 25 /im, 
and the infrared background expected at this wavelength. 

Soifer & Neugebauer (1991) calculated the local luminosity density in the four IRAS 
passbands based on their derived luminosity functions. At 25 /im, they found that the log 
of the luminosity density in units of L Q Mpc -3 is 7.2. Integrating our parametric function 
gives a log luminosity density of 6.9, which is a factor of 2 lower in linear terms than 
the Soifer & Neugebauer result. The differences in the sub-L* portion of the luminosity 
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functions account for the different integrated luminosity densities. 

Hacking & Soifer (1991) presented number count distributions at 25 /im, using a fit to 
the Soifer & Neugebauer (1991) luminosity function. (Those authors fit a hyperbola to 
the visibility function, which has the power-law asymptotes at high and low luminosities 
like our parametric form, except that it forces the absolute values of the slopes of the 
asymptotes to be identical.) As they noted, the Soifer & Neugebauer functions at other 
wavelengths are 20%-30% higher than those derived from other samples (Saunders et al. 
1990). Our improved 25 fim function is similarly about 15% lower in normalization than 
the Soifer & Neugebauer estimate. Hence, models based on this luminosity function and its 
normalization will produce number counts about 15% lower at the bright end than those in 
Figure 1(c) of Hacking & Soifer (1991). Such a change would bring the number counts in 
better agreement with the IRAS number counts presented in their same figure. 

Finally, we have substituted our luminosity function into the backgrounds model of 
Hacking & Soifer (1991). The modified model yields estimates of the infrared background 
at 25 jum that are factors of 2.5 to 2.7 times smaller than the original Hacking & Soifer 
estimates. Our results indicate that the detection of an IR background at this wavelength 
is likely to be even more difficult than previously predicted. 



6.4. Evolution 

To this point, we have neglected any effects of galaxy evolution on our derivation of the 
luminosity function. Since the maximum-likelihood estimators are independent of density 
variations, the shape of the luminosity function will be unchanged by density evolution 
of the galaxy population. Our normalization considers galaxies out to 20,000 km s™ 1 , a 
redshift where evolution effects should be small. 
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More importantly, our maximum likelihood fit to the radial density (section |5| and 
Figure ^) does not show a systematic increase with redshift and hence no indication of 
density evolution. This result is in direct contrast to that of Saunders et al. (1990). From 
their 60-/iin-selected sample, those authors found that density evolution of (1 + z) e ' 7 
matched their radial density fit. However, we note that the existence of evolution from 
IRAS samples is not yet settled — for the deepest IRAS sample, Ashby et al. (1996) did 
not find the high-redshift tail expected from strong evolution. 

Nevertheless, we have computed the shape of the parametric luminosity function for 
the case of exponential luminosity evolution corresponding to (1 + z) 3 at low redshifts (cf. 
Saunders et al. (1990)), by correcting the luminosity of each galaxy to the present epoch. 
The values of the shape parameters of this function are a = 0.437±0.032, (3 = 1.769±0.067, 
and = 3.98 ± 0.34 x 10 9 L o . The only change from the function described in Table |3] is 
that the slope of the luminosity function at the high-luminosity end is slightly lowered. 



7. Summary and conclusions 

The following are the results of this paper: 

1. We have selected a sample of 1458 galaxies with redshifts from the IRAS Faint 
Source Survey with a flux density limit of 250 mJy at 25 fim. An additional 17 galaxies do 
not have redshifts available. 

2. The local luminosity function is derived using the 1/V max estimator and both 
parametric and non-parametric maximum likelihood methods. The 1/V max estimate is 
significantly affected by the Local Supercluster. The maximum likelihood methods are 
independent of density variations, and we consider the parametric fit with parameters in 
Table |3| to be the best estimate of the local luminosity function at 25 /im. 



- 21 - 



3. A maximum likelihood fit to the radial density in this sample is used to correct the 
^/Vmax estimate. The fit shows no sign of a systematic increase with redshift of the density, 
as would result from density evolution of the galaxy population. 

4. The 1/V max luminosity function derived from a smaller sample by Soifer & 
Neugebauer (1991) is significantly contaminated by the Local Supercluster. Predictions of 
number counts and local luminosity density based on that function are 15-20% higher than 
those indicated by our improved luminosity function. The new function also leads to lower 
predictions of the mid-infrared background due to galaxies. 

We thank Susan Tokarz, Perry Berlind, and Jim Peters for observations made at F.L. 
Whipple Observatory. We also thank Iffat Khan and Joe Mazzarella for assistance with 
database searches, and Carol Lonsdale and Tom Soifer for helpful discussions. We are 
grateful to Cong Xu for a critical reading of this manuscript. We thank the anonymous 
referee and our Scientific Editor, Steven Willner, for numerous comments that have 
improved this paper. The SIMBAD database is maintained by CDS in Strasbourg, France. 
The NED database is supported at IPAC by NASA. 

A. The Distance to Virgo and a Different Hubble Constant 

We have assumed a distance 17.6 Mpc to Virgo in the linear Virgocentric inflow model 
(Aaronson et al. 1982) to calculate the distance to each galaxy which does not have a 
secondary distance. To be consistent with our assumption of 75 km s _1 Mpc -1 for the 
Hubble constant, we have used 1019 km s _1 as the observed velocity of Virgo corrected 
to the centroid of the Local Group, and implicitly assumed an infall velocity of about 300 
km s _1 for the Local Group toward Virgo. The resulting Hubble flow velocity of Virgo, 
~1320 km s _1 , is consistent with many other estimates (e.g. Faber et al. 1989; Tonry et 
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al. 1990; Mould et al. 1995; Huchra 1996). A different Hubble constant would change the 
distance estimate to Virgo if the Hubble flow velocity of Virgo remains unchanged. 

We have calculated the l/V max luminosity function assuming no Virgocentric inflows, in 
which we used the redshift as the direct distance measure unless a secondary distance exists 



for a galaxy. Figure [14| compares the 1/V max luminosity functions for the 250 mJy-limited 
sample with and without Virgocentric inflow. A Hubble constant of 75 km s _1 Mpc -1 is 
used in both cases. The two estimates agree very well. This indicates that the inflow model 
we use has a small impact on the estimate of the distance of the galaxies in the sample, 
and that most of these galaxies are not strongly affected by the gravity from Virgo. Note 
that a luminosity function simply scales with the Hubble constant (the luminosity changes 
as Hq 2 ) if all distance estimates in a sample are from redshifts. 

The number of galaxies with secondary distances in our sample is small (only 22). 
We have done experiments to verify that these galaxies do not contribute statistically 
significantly to the calculated luminosity function. We therefore conclude that our results 
can be simply scaled if a different Hubble constant is used. 



-23- 
REFERENCES 

Ashby, M.L.N., Hacking, P.B., Houck, J.R., Soifer, B.T., & Weisstein, E.W. 1996, ApJ, 465, 
428. 

Aaronson, M., Huchra, J., Mould, J., Schechter, P. L., and Tully, R. B. 1982, ApJ, 258, 64. 
Boulade, O. et al. 1996, A&A, 315, L85. 

Chokshi, A., Lonsdale, C. J., Mazzei, P., and De Zotti 1994, ApJ, 424, 578. 
Choloniewski, J. 1986, MNRAS, 223, 1. 
Davis, M., & Huchra, J. 1982, ApJ, 254, 437. 

Efstathiou, G., Ellis, R. S., & Peterson, B.A. 1988, MNRAS, 232, 431. 

Egret, D., Didelon, P., McLean, B. J., Russell, J. L., and Turon, C. 1992, A&A, 258, 217 

Faber, S. M., Wegner, G., Burstein, D., Davies, R. L., Dressier, A., Lynden-Bell, D., and 
Terlevich, R. J. 1989, ApJS, 69, 763. 

Fang, F., Shupe, D.L., Xu, C, and Hacking, P.B. 1998, ApJ, submitted. 

Fairall, A., Lowe, L. & Dobie, P. 1988, Pub. of the Department of Astronomy of the 
University of Cape Town, No. 5, A Catalog of Galaxies South of Declination -30 
That Have Been Observed Spectroscopically. 

Felton, J. E. 1977, AJ, 82, 861. 

Fisher, K. B., Huchra, J.P, Strauss, M., Davis, M., Yahil, A., & Schlegel, D. 1995, ApJS, 
100, 69. 

Hacking, P. B. et al. 1996, The Wide-Field Infrared Explorer (WIRE) Mission, to appear 
in International Symposium on Diffuse Infrared Radiation and the IRTS, Institute 
of Space and Astronomical Sciences, Sagamihara, Japan. 

Hacking, P. B., Condon, J. J., and Houck, J. R. 1987, ApJ, 316, L15. 



-24- 



Hacking, P. B., and Soifer, B. T. 1991, ApJ, 367, L49. 

Huchra, J., Davis, M., Tonry, J. & Latham, D. 1983, ApJS, 52, 89. 

Huchra, J. P. 1996, in IAU Symposium No. 168, Galaxy Motions in the Nearby Universe, 
pp. 143. 

Kendall, M.G., Stuart, A., & Ord, J.K. 1987, The Advanced Theory of Statistics, Vol. 2 
(5th ed.; New York; Oxford University Press). 

Knapp, G. R., Rupen, M. P., Fich, M., Harper, D. A., Wynn- Williams, C. G. 1996, A&A, 
315, L75. 

Lasker, B. M., Sturch, C. R., McLean, B. J., Russell, J. L., Jenkner, H., & Shara, M. M. 
1990, ApJ, 99, 2019 

Lawrence, A., Rowan-Robinson, M., Crawford, J., Parry, L, Xia, X-Y., Ellis, R.S., 

Frenck, C.S., Saunders, W., Efstathiou, G. and Kaiser, N., 1994 MNRAS to be 
submitted;preprint 

Longair, M. S. 1978, Observational Cosmology, A. Maeder, L. Martinet, and G. Tammann 
(Geneva Observatory). 

Lynden-Bell, D. 1971, MNRAS, 155, 95. 

Marzke, R.O., Geller, M.J., Huchra, J. P., & Corwin, H.G. 1994, AJ, 108, 437. 

Moshir M., et al. 1992, IRAS Faint Source Catalog Explanatory Supplement, Version 2. 

Mould, J. et al. 1995, ApJ, 449, 413. 

Rowan-Robinson, M. et al. 1996, Spectral Energy Distributions, Starburst Models and Star 
Formation History, MNRAS, submitted. 

Sandage, A., Tamman, G. A., & Yahil, A. 1979, ApJ, 232, 352. 

Saunders, W., Rowan-Robinson, M., Lawrence, A., Efstathiou, G., Kaiser, N., Ellis, R. S., 
& Frenk, C. S. 1990, MNRAS, 242, 318. 



-25 - 

Schember, H., Kemp, J., Ames, H., Hacking, P., Herter, T., Fafaul, B. Everett, D., Sparr, 
L. 1996, "The Wide-Field Infrared Explorer (WIRE)", SPIE Proceedings, 2744, 
Infrared Technology and Applications XXII, p. 751-760. 

Schmidt, M. 1968, ApJ, 151, 393. 

Soifer, B. T. et al. 1987, ApJ, 320, 238. 

Soifer, B. T., & Neugebauer, G., 1991, AJ, 101, 354. 

Strauss, M., Davis, M., Yahil, A., & Huchra, J. P. 1990, ApJ, 361, 49. 

Tonry, J., Ajhar, E., & Luppino, G. 1990, AJ, 100, 1416. 

Yahil, A., Strauss, M. A., Davis, M., & Huchra, J. P. 1991, ApJ, 372, 380. 

Yentis et al. 1992, in "Digitised Optical Sky Surveys", p. 67, ed. MacGillivray & Thomson 
(Dordrecht : Kluwer ) . 



This manuscript was prepared with the AAS IATjtX macros v4.0. 



-26- 



Fig. 1. — Color-color plots for the galaxies with redshifts in the sample. The upper plot 
shows galaxies with good or moderate detections in each of the 12, 25 and 60 /im bands. 
The lower plot show galaxies with upper limits in either or both the 12 or 60 /im bands. The 
hatched region shows the region of the color-color diagram that is excluded from the sample 
selection. 

Fig. 2. — Differential number counts versus flux for the galaxies with redshifts in the sample. 
The line has a slope of -1.5. 

Fig. 3. — Visibility plots for the 250 mJy sample. The open triangles represent the results of 
the 1/V max estimate of the luminosity function fit to the entire sample. The solid line is the 
parametric maximum likelihood estimate for galaxies with redshifts greater than 500 km s _1 . 
The parameters of the luminosity function are listed in Table [| The curve is truncated at 
L s = 4 x 10 7 L Q . The distance scale at top is computed from the 250 mJy flux density limit 
and the luminosity scale at the bottom. 

Fig. 4. — V/V max for the full sample (upper panel), and separately for northern and southern 
Galactic latitudes (lower panel). The dashed line at V/V max = 0.5 shows the expected value 
for a uniform distribution of galaxies. 

Fig. 5.— The non-parametric maximum-likelihood estimate of the visibility function 
(circles). The parametric maximum-likelihood estimate from Figure ^| is also shown for 
comparison. The different estimators agree within the la error bars shown on the Figure. 

Fig. 6. — The l/V max visibility functions for the high-quality (open squares) and low-quality 
subsamples (filled squares) defined in Section The la error bars are shown. 

Fig. 7. — The 1/V max visibility functions for the high-quality subsample (filled-squares) 
defined in Section |]and the full 250 mJy-limited sample (filled-circles) . Also shown are the 
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la error bars. 

Fig. 8. — Distribution of the quantity ^/(L)/\l/(L m j n ) in the sample for galaxies in different 
distance ranges. The mean value is normalized to 1 in each distance range. The flatness of 
these distributions confirms that the completeness of the sample is high to at least 400 Mpc. 

Fig. 9.— The non-parametric maximum-likelihood estimate (solid histogram) of the 
radial density distribution for the 250 mJy flux-limited sample at 25 /im. This estimate 
is independent of the galaxy luminosity function. The dashed lines show the observed radial 
density distribution. 

Fig. 10. — The V/V max plot for the full flux-limited sample after correcting for density 
variations as described in the text. The la error bars are shown. 

Fig. 11. — The density-corrected 1/V max estimate of the visibility function. The parametric 
maximum-likelihood estimate from Figure § is also shown for comparison. The density 
inhomogeneities have been corrected according to the radial density estimate shown in Figure 
|9|. The different estimators agree within the la error bars shown on the Figure. 

Fig. 12. — The observed redshift distribution (dashed-line histogram) and the distribution 
corrected for the density variations in Figure |9] (solid-line histogram) as described in the 
text. The dotted line is the distribution predicted by the parametric maximum-likelihood 
luminosity function from Figure [5] and Table 

Fig. 13.— Local luminosity functions at 25 /im plotted in terms of space density as a 
function of luminosity. The dashed lines were computed by varying the values in Table |3| by 
their la uncertainties. 

Fig. 14. — The 1/V max visibility functions for the 250 mJy-limited sample with (filled circles, 
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as in Figure [7]) and without (circles) the Virgo-centric inflow. The two functions agree well 
within the la error bars, indicating that most of the galaxies in the sample are not affected 
by the inflow model. 
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Table 1. Source of redshifts for sample 



Source of redshift Number 



Public Sources: 

1.2 Jy catalog 1277 

Nov. 1993 ZCAT 49 

NED 79 

SIMBAD 6 

Fairall, Lowe, & Dobie (1988) 1 

Private Sources: 

FLWO 1.5m (Table |) 41 

QDOT (Lawrence et al. 1994) 5 
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Table 2. Previously unpublished redshifts from FLWO 



R.A. Dec. F„(25) cz 

(B1950) (Jy) (km s" 1 ) 



00 


17 


02 


.7 


-04 


56 


46 





.3734 


6181 + 28 


00 


20 


47 


.6 


+10 


29 


50 


0. 


.2633 


69094 ± 85 


00 


37 


17. 


.8 


+08 


29 


18 


0. 


.2619 


17482 ± 34 


00 


44 


30 


.0 


-18 


03 


33 





.2776 


43979 ± 48 


00 


44 


37 


.2 


+ 10 


14 


46 





.3265 


49959 ± 54(iV) 
50133 + 43(5) 


00 


56 


04 


.2 


+08 


31 


54 





.2861 


17269 ± 36 


01 


27 


44 


.3 


+07 


53 


04 





.2541 


33899 ± 44 


01 


47 


43. 


.0 


-16 


55 


12 





.2875 


48432 ± 39 


07 


55 


56. 


.0 


+50 


58 


36 





.3589 


16305 + 41 


08 


11 


42. 


.0 


+46 


13 


00 





.2951 


12277+ 136 


08 


46 


21 


.1 


+ 11 


26 


03 





.2626 


23218 + 33 


08 


57 


06 


.0 


+18 


01 


00 





.2645 


15940 ± 20(E) 
15830 + 19(W) 


09 


39 


16 


.1 


+23 


54 


38 





.3299 


6407 ± 39 


09 


42 


35 


.4 


+ 17 


51 


45 





.4363 


38440 ± 41 


10 


03 


12. 


.0 


+13 


12 


00 





.4967 


2789 ± 28 


10 


07 


43. 


.0 


+14 


27 


39 





.2729 


59689 ± 55 


10 


08 


05. 


.4 


+06 


26 


45 





.2883 


29445 ± 39 


10 


23 


48. 


.0 


+00 


56 


00 





.3076 


6422 ± 18 



Table 2 — Continued 



R.A. Dec. F„(25) cz 

(B1950) (Jy) (km s" 1 ) 



11 


06 


26. 


.0 


-23 


21 


57 





.3648 


24551 ± 41 


11 


19 


58. 


.7 


+04 


31 


06 





.5010 


11399 + 38 


11 


22 


24. 


,9 


+02 


56 


02 





.2533 


18813+ 17 (N) 
18857 + 24(5) 


11 


31 


56. 


.3 


+22 


47 


49 





.3011 


9230 ± 34 


11 


57 


29. 


.2 


-03 


30 


20 





.2663 


8024 ± 20 


12 


12 


17. 


.0 


-03 


12 


00 





.4950 


10108 + 37 


12 


17 


18. 


.6 


+09 


12 


59 





.2743 


7398 ± 31 


12 


19 


41. 


.7 


+00 


09 


03 





.2991 


51652 ±29(JV) 
olozo ± ob{o ) 


12 


23 


25. 





-23 


19 


36 





.3232 


14568 ± 40 


12 


36 


16. 


.3 


-27 


02 


00 


0. 


.2678 


7502 ± 22 


12 


38 


18. 


.7 


+27 


50 


19 





.3110 


16947 ± 34 


12 


50 


30 


.4 


-02 


56 


36 


0. 


.4896 


6807 ± 34 


13 


24 


15. 


.4 


-06 


08 


38 





.2963 


13489 ± 46 


13 


53 


28. 


3 


-07 


51 


36 





.3494 


22750 ± 40 


13 


55 


57. 


.1 


-15 


53 


51 


0. 


.3377 


10868 ± 34 


15 


19 


54. 


4 


-13 


51 


29 





.2555 


7182 ± 25 


15 


21 


47. 


.0 


-11 


02 


07 





.2667 


15447 ± 37 


15 


21 


48. 


.0 


+05 


01 


00 





.2588 


10715 + 10 
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Table 2 — Continued 



R.A. 


Dec. 


F v (25) 


cz 


(B1950) 


(Jy) 


(km s -1 ) 


21 11 45.0 


-27 31 29 


0.2655 


25801 ± 33 


22 34 59.5 


-25 44 16 


0.2646 


23353 ± 21 


23 13 02.6 


-00 58 50 


0.3414 


8318 ±36 


23 29 00.6 


-12 52 29 


0.3778 


6357 ± 40 








11478 ±23(W) 


23 40 45.2 


+04 16 37 


0.3489 


11564 ±21(iV) 








11334 ±21(5) 



Note. - Multiple redshift components are labeled 
with letters denoting their relative orientations. 
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Table 3. Fitted Parameters of Luminosity Function 



a 






C 








(Le) 


(lir 3 M pc - 3 ) 


(10- 3 Mpc" 3 ) 



0.437 ±0.032 1.749 ±0.067 4.07 ± 0.36 x 10 9 



1.23 ±0.09 



9.2 ±0.7 
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